Doppler-free Adiabatic Self-Induced Transparency 
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We demonstrate that a Doppler broadened two-level medium can be made transparent to a laser 
pulse by an appropriate adiabatic variation of the laser field amplitude and its nominal detuning. 
This new technique of adiabatic self-induced transparency (ASIT) is compared with the well known 
self-induced transparency (SIT) phenomenon, showing that the adiabatic method is much more ro- 
bust to variations of the system parameters. We also discuss a possible experimental implementation 
of ASIT using 87 Rb atoms. 
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I. INTRODUCTION 

The ability to tailor the optical properties of a medium 
in such a way that it becomes transparent to a reso- 
nant laser field while its index of refraction is shaped at 
wish has been used for a wide range of applications, e.g., 
from slow light [l], Q and dark-state polariton physics Q 
to light-matter interfaces These investigations are 
mainly motivated by the fast development of quantum 
information technologies and have been carried out in 
different types of media from atomic vapours, either cold 
1 1 or at room temperature Q , to solid state materials 
E 11] > photonic crystals @, [I[ and optical fibers @ . 

For practical implementations it would be desirable to 
work at room temperature but often this implies the pres- 
ence of inhomogeneous broadening that tends to reduce 
the transparency. Explicitly, in hot atomic vapors, the 
atomic motion leads to Doppler broadening due to the 
velocity-dependent shift of the laser field frequencies. If 
so, one direct approach to avoid inhomogeneous broad- 
ening consists of using cold atomic samples although it 
involves demanding technological requirements. In fact, 
albeit there are some proposals on how to reduce the neg- 
ative role of the Doppler broadening in a hot vapor [lCJ , 
a fully satisfactory solution is still challenging. Here we 
will focus on avoiding the adverse effects of the Doppler 
broadening to obtain on-resonance transparency for two- 
level hot atomic vapours. However, our results could be 
extended to other inhomogeneously broadened systems, 
e.g., solids with non-uniform crystalline fields [111 ]. 

In a two-level system, there exists a simple approach 
for achieving Doppler-free transparency named self- 
induced transparency (SIT) [l~2j], which consists in shap- 
ing a resonant optical pulse propagating in the medium 
such that an integer number of Rabi oscillations is per- 
formed. This occurs in the coherent regime whenever the 
pulse area is an integer multiple of 2ir. Therefore, SIT 
requires a very precise temporal control of the Rabi fre- 
quency. Analogously to SIT, when the pulse area is a 
multiple of an odd number of 7T, self-induced absorption 



(SIA) occurs. 

An alternative approach to achieve resonant trans- 
parency consists in using induced atomic coherences in 
three- or multi-level atomic systems, e.g. , by means of the 
coherent population trapping (CPT) [13[ and the electro- 
magnetically induced transparency (EIT) [l4| phenom- 
ena or via a double-STIRAP process [Hj]. In most cases, 
an intense laser field, the so-called control field, couples 
to an atomic transition while modifying the optical prop- 
erties of the medium for a weak laser field, the probe, 
acting on an adjacent transition. However, the fact that 
all these atomic coherence effects, CPT, EIT, STIRAP 
and doublc-STIRAP, operate in the so-called two-photon 
(or Raman) resonance condition means that, in the pres- 
ence of Doppler broadening, the two laser fields must 
have near identical frequencies which severely restricts 
the range of applications of atomic coherence effects in 
hot atomic vapors (loj . 

In this paper, we present a novel approach that gives 
rise to Doppler-free transparency in a two-level (2L) sys- 
tem resembling a doublc-STIRAP process in a three-level 
(3L) system. To achieve self-transparency we will appro- 
priately modify in time the amplitude and detuning of the 
laser pulse. It will be shown that this technique, which 
we name adiabatic self-induced transparency (ASIT), op- 
erates for a larger parameter range than SIT and is much 
more robust against fluctuations of the parameter values. 

The paper is organized as follows: in Section [H] the na- 
ture of the ASIT phenomenon is described in the frame- 
work of the topological equivalence between density ma- 
trix equations for a 2L system and the Schrodinger equa- 
tion for a 3L system. In Section IIII1 the time evolution 
of the 2L system under conditions for ASIT is studied in 
the inhomogeneously (Doppler) broadened case. Section 
HVl is devoted to the comparison between the proposed 
ASIT process and the already known SIT phenomenon. 
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II. ADIABATIC SELF INDUCED 
TRANSPARENCY 

We consider a two- level atom (Fig. Q] (a)) interact- 
ing with a classical electromagnetic field (EM), e(t) = 
eE(t) cos(ui) in the electric dipole approximation. Under 
the rotating wave and slowly varying envelope approxi- 
mations, the coherent dynamics of the two-level atom is 
given by the following density matrix or optical Bloch 
equations: 



dt \ w 



A 
-A -CI 
CI 




(1) 



where we have introduced the real-valued variables that 
define the coordinates of the Bloch vector moving in a 
3D space 

U = 2Rc(a 12 ), y = 2Im(ai 2 ), W = (a 22 - a n ) , 

an and cr^ with i, j = 1, 2 account for the corresponding 
populations and coherences. A = lj— u>2\ is the frequency 
detuning of the field with respect to the instantaneous 
atomic transition frequency w 21 . fl(t) = — dy2.eE (t) /h 
denotes the time dependent Rabi frequency with c? 12 be- 
ing the electric dipole moment of the atomic transition. 
The Rabi frequency, without loss of generality, is taken 
real and positive. Throughout the paper, we will consider 
that both the detuning and the optical field envelopes are 
time dependent. We name the temporally shaped detun- 
ing as D pulse, and the temporally shaped amplitude field 
as EM pulse. 

Vitanov et al. [T3| have recently shown that Eqs. ^ 
are topologically equivalent to the Schrodinger equation 
for a three-level atom in a A configuration interacting 
with two resonant laser fields (Fig. Q] (b)) under the fol- 
lowing identifications: U <^> C3, V <3>_C%, W 4^ C\, 
and 2 A (t) <^> fl s (£), 2Q Q P , where Ci,3 = Ci, 3 and 
C-2 = — iC 2 , with Ci being the probability amplitude to 
find a three-level atom in state Clp^s denote the Rabi 
frequencies of the two applied fields. One should remark 
that for a 2L system the variables U, V and W are real 
quantities, while for a 3L system the variables Ci are, in 
general, complex numbers. The analogy between the 2L 
and the 3L systems is shown in Fig. [ljc). 

In a 3L system, STIRAP is a widely used technique for 
transferring the atomic/molecular population between 
the two lower levels of a A-scheme and it is based on 
adiabatically following one of the three energy eigen- 
states of the system given by |3>(f)) = cos9(t)|I) — 
sinO(i)|3), with the mixing angle defined as Q(t) = 
arctan [Tip (t)/ris(t)]. In a 2L system and according to 
the previous analogy, the equivalent to the former en- 
ergy eigenstate is the following combination of the pop- 
ulation difference and the real part of the atomic coher- 
ence 4>{t) = cos9{t) x W(t) - sm9(t) x U(t), with the 
mixing angle being now 0(t) — arctan [Cl(t)/A(t)]. In- 
deed, this analogy opens the possibility to directly extend 
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FIG. 1: (Color online) (a)Two-level system under considera- 
tion where Q and A are the Rabi frequency and the field-atom 
detuning; (b) Three-level system interacting with two fields 
with Rabi frequencies fip and Qs', (c) equivalence between 
the density matrix variables of the two-level system and the 
three-level picture. Schematic representation, for the two- 
level case, of the Bloch vector evolution on the Bloch sphere 
of the (U, V, W) variables, under (d) SIT with 2tt pulses, (e) 
STIRAP and (f) ASIT processes. 



the results that were known previously for 3L systems, 
obtained from the adiabatic temporal variation of the 
Rabi frequencies of the two coupled fields, to 2L systems, 
by appropriately engineering the temporal profile of the 
detuning and the Rabi frequency of a single EM pulse. 
We note that the temporal variation of the field ampli- 
tude and the detuning can be accomplished for instance 
by chirped pulses, and that chirped pulse excitation to 
achieve population inversion in a 2L system has been 
discussed via Stark Chirped Raman Adiabatic Passage 
(SCRAP) [lj. 

Such analogy between STIRAP and chirped two-state 
excitation has been studied in (l7j in the complete co- 
herent case, e.g., without spontaneous emission, and in 
the absence of Dopplcr broadening. From that study it 
follows that for the 2L-STIRAP case the Bloch vector 
moves from the south pole (W = — 1, U = and V = 0) 
to the equator (W = 0, U = 1 and V = 0) [see Fig. [He)] 
on the Bloch sphere. Starting with the population in 
the lower level and by applying first a D pulse and later 
on, and with an appropriate overlap, the EM pulse, the 
system ends up in a coherent superposition of the two 
states. In this process, the variable V accounts for EM 
field absorption (positive values) and stimulated emission 
(negative values), U describes the dispersion properties 
of the transition, and W gives the transfer of population. 

In a 2L system, a functional definition of trans- 
parency is that the two-level atom evolves coherently 
from Wifinitiai) = ~ 1 (population in the lower level) 
back to W(t fi na i) = — 1. This coherent evolution should 
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be fast enough such that spontaneous emission becomes 
negligible, which in turn, means that there is not irre- 
versible EM absorption. Such definition of transparency 
is equivalent to state that the time averaged imaginary 
part of the atomic coherence weighted by the Rabi fre- 
quency is equal to zero, i.e., g (ftU) = 0, where g = 

N 2hloc 2 * s ^ ne S am parameter per unit length with N 
the atomic density, eo the electric permitivity and c the 
speed of light. Hereafter the brackets (...) mean time 
averaging. 

As it has been mentioned above, a well known method 
to obtain transparency in a 2L system is the SIT phe- 
nomenon (l2l |. i.e., coherent induced absorption of pulse 
energy during the first half of the pulse followed by co- 
herent induced emission of the same amount of energy 
back into the beam direction during the second half of 
the pulse. In this case, the complete transfer of the 
population from the lower to the upper state and back 
can be realized with a fully resonant excitation whenever 
the pulse area is a multiple of 27r. Pulse areas multi- 
ple of an odd number of ir result in complete absorption 
(SIA). The SIT/SIA phenomena have an absorptive na- 
ture, which means that in the Bloch sphere the evolution 
occurs predominantly in the plane (V, W) with {7 = 
[see Fig. H^d)], i.e. the transfer of the population is due 
to the EM absorption and emission processes involved 
and, therefore, g (QV) = with g 2 ((ftU) 2 ) / (ft 2 ) 2 ^ 0. 

The novel method to obtain Doppler-free transparency 
in a 2L system, which we propose here, i.e. ASIT, con- 
sists in a process that resembles a double-STIRAP se- 
quence in a 3L system. In a 3L system (Fig. [Tfb)), a 
double-STIRAP process starts and ends in the ground 
state |1) since the mixing angle starts at 6 = 0, goes 
to 6 = 7r/2, and ends back again at = 0. To pro- 
duce a double-STIRAP process, a first Stokes pulse has 
to be applied at the beginning of the process, then a 
partially overlapping pump pulse, and finally a partially 
overlapping second Stokes pulse. Alternatively, a rela- 
tively longer Stokes pulse that overcomes the pump pulse 
can be also used. In a similar manner, ASIT can be pro- 
duced in a 2L system by applying an EM pulse in be- 
tween two D pulses. The ASIT phenomenon is mainly 
of dispersive nature, since the evolution occurs predom- 
inantly in the plane (U, W) , with much weaker relative 
excitation of the variable V related with absorption and 
of variable W [see Fig. Q] (f)] and, therefore, g (QV) = 
and g 2 ((ftU) 2 ) / (ft 2 )" takes values smaller than in the 
SIT case. 

Up to now we have not considered the incoherent pro- 
cesses that might be present in the system. In fact, there 
are no analogies between the 2L and 3L systems regard- 
ing incoherent processes like spontaneous emission or de- 
phasing. In the 3L case, the spontaneous decay rate of 
the excited state does not play any role either in the 
STIRAP or in the double-STIRAP processes, because 
the whole adiabatic process occurs without coupling the 
excited state. In contrast, in the 2L case, spontaneous 
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FIG. 2: (Color online) Temporal profiles of the EM (black 
line) and D (magenta line) pulses (a,c,e,g). Time evolution 
of the population difference (W) and the medium coherence 
(U, V) (b,d,f,h) in a homogeneously broadened two-level sys- 
tem under the SIA (a,b), SIT (c,d) and ASIT (e-h) processes. 
The areas of the pulses are equal to A = An = Aa = 9-7T 
(a,b,e,f) and = 10-7T (c,d,g,h). The D pulses and the rise and 
fall parts of the EM pulses have gaussian profile. The tem- 
poral width of the D pulses is ta = 0.0037 -1 , Ao = 0. 
A/7 = 5318 (e,f) and = 5909 (g,h). tl /y = 1772 (a,b,e,f) 
and = 2363 (c,d,g,h). The rise 7> and fall 77 time of the EM 
pulses is equal to r r j — ta/2. The D pulses are shifted away 
by At = 1.5ta from the top flat part of the EM pulse. The 
time duration of the top flat part of the EM pulses is set to 
keep the EM pulse area indicated above. In all cases repre- 
sented, the detuning of the EM pulses during the top flat part 
is close to zero value. 



emission plays an important role because the exited state 
is always involved in the process. We will assume that 
the population difference decays at a rate 7 while the de- 
cay of the real and imaginary parts of the coherence is 
r/2. In general, T accounts for several incoherent pro- 
cesses like spontaneous emission, collisions, dephasing or 
phase fluctuations. Here we only consider the effect of 
the spontaneous emission in the decay of the coherences 
and therefore T = 2j. The requirement of full coher- 
ence for ASIT to work imposes a maximum time duration 
of the process limited by the lifetime of the levels, i.e., 
T7 <C l, where T is the total time for the ASIT process. 
Under these conditions, comparison between the three 
just above described processes of SIA, SIT and ASIT is 
shown in Fig. [2] where the time evolution of the variables 
(W, U, V) in a homogeneously broadened 2L medium is 
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represented. It should be mentioned that in the case we 
consider here Figs. [2] (b,d,f,h) represent the temporal re- 
sponse of a single 2L atom as well as of a homogeneously 
broadened medium composed of such 2L atoms. In con- 
trary to the SIA and SIT processes (see Figs. [2] (a-d)) 
where the final state of the system strongly depends on 
the EM pulse area, under the ASIT phenomenon the pop- 
ulation returns back to the lower level, giving rise to the 
transparency of the 2L medium, independently on the 
area of the exciting EM pulse as demonstrated in Figs. [5] 
(c-h). Although the strength of the light-atom interac- 
tion is larger for the SIA and SIT processes than for the 
ASIT process, but it is strong in all cases. 



III. ASIT IN INHOMOGENEOUSLY 
BROADENED MEDIA 




Maxwell's 
distribution 
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Now we focus on the role of the Doppler broaden- 
ing in the 2L system under consideration by assuming a 
Maxwell distribution for the projection of the atomic ve- 
locities on the light propagation direction. We will solve 
the ASIT process for different atomic velocity classes and 
sum up the total response along the Maxwell distribu- 
tion. For an atom moving with some velocity v along the 
direction of the propagating laser field (with wavenum- 
bcr k = u/c), the instantaneous field detuning will be 
A (v, t) = A (t) — kv , where A (t) is the nominal detun- 
ing for atoms at rest. For dcfmitcness we choose temporal 
Gaussian profiles for the D and the EM pulses: 

= n e- t2/T n , (2) 

where A defines the amplitude of the temporal variation 
of the detuning with respect to the constant detuning 
offset A , r2 is the amplitude of the EM pulse, ta,q are 
the time durations of the D and EM pulses, and ±Ar 
are the delays of the D pulses with respect to the center 
of the EM pulse. 

For a complete transfer of the atoms, the double- 
STIRAP process should be performed adiabatically, con- 
dition that for atoms moving with velocity v reads [13]: 

T^j fijj + A 2 (v,t) > 10 . (3) 

Let us assume, for simplicity and without losing gen- 
erality, that Ao and A are both positive. Then for the 
homogeneous case or for a small Doppler width, i.e. for 
Aojd < A, the adiabaticity condition ([3]) can be fulfilled 
for all atoms simultaneously, provided that appropriate 
values for the parameters, in particular for Oo, are cho- 
sen. This might not be the case, however, for hot atomic 
vapours, where the velocity-dependent Doppler shift can 
strongly modify the D pulse offset during the ASIT pro- 
cess for some atomic velocity classes. The largest reduc- 
tion in the instantaneous values of the detuning A 2 (v, t) 



FIG. 3: (Color online) Distributions for the final value (red 
solid line) of the atomic population inversion, Wfi n , and of 
the maximum (empty circles) and minimum (crosses) values 
for the mixing angle 9 with respect to the Doppler shift kv 
at different amplitudes of the EM and D pulses (b-d). The 
areas of the EM and D pulses, defined as An = QoTn/V2 and 
A a = Ata/\/2, are equal to each other A — Aq/tv — Aa/k 
= 2 (b), = 6 (c), = 10 (d). The Maxwell distribution of the 
atomic velocities is shown in (a) (dashed line). The inset in 
(a) presents the time evolution of the exciting EM (fi) and D 
(A) pulses. Abrupt decrease in the 6 max corresponds to the 
velocity class for which kv ~ Ao + A. The temporal widths of 
the input pulses are ta = rn = O.OI7 -1 , Ao = and At = 
1.5rn. The Doppler width is fixed at Alod = 10 3 7 and V = 
27. Macroscopic response of the inhomogeneously broadened 
2L medium is given by the integration of the atomic response 
[shown by red solid lines in figures (b-d)] weighted by the 
probability density for the atomic velocities of the Maxwell 
distribution [shown by the dashed line in figure (a)]. 



will occur for atomic velocities verifying kv > Aq and in 
particular for kv ~ (Ao + A), where the detuning will re- 
main close to zero during significant time intervals (those 
where the D pulse passes near its maxima values). This 
can break the adiabaticity condition ([3]) for atoms with 
such velocity, so that they can be transferred to the ex- 
cited state breaking the transparency under the doublc- 
STIRAP process. One needs to note that to maintain 
adiabaticity of the process for small detuning A (v, t) the 
Rabi frequency Qp should be increased substantially, see 
for instance [l3 ]. 

We have studied in detail the response of 2L atoms 
with different values of the Doppler shift kv along the 
Maxwell distribution by exciting them with D and EM 
pulses of the form given in . Some of the results are 
presented in Figs. [3] and 2J for a specific Doppler fre- 
quency profile and a specific temporal pulse profile [see 
Fig. [3J (a)]. Fig. 5] shows the total atom-field detuning 
(sum of the D-pulsc value plus the Doppler shift kv) [see, 
Figs. |U (a), (c) and (e)], and the atomic-state evolution 
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FIG. 4: (Color online) Temporal profiles for the EM pulse 
(black line) and the total detuning seen by the atom (magenta 
line), for three different atomic velocity classes (a,c,e): kv — 
Awd3/4 (a,b), kv — Awd/3 (c,d) and kv = — Aojd/3 (e,f). 
Figs, (b), (d) and (f) represent the time evolution of the 
atomic variables W, U and V. The pulse areas are An ~ 
Aa = Wit. The other parameters are the same as in Fig. [3] 



for three different atomic velocity classes v [see, Figs. 0] 
(b), (d) and (f)]. It is clear from Fig. [4] that in all cases 
the light-atom interaction is quite strong [it is less strong 
in the case of Figs. [4] (e) and (f ) because in that case the 
total field detuning is largest] , but always the final state 
of the population difference W is the same as the initial 
state, i.e., W = — 1. This is due to the adiabaticity of 
the ASIT process. The dependence of the final value of 
the atomic population inversion, Wfi n , i.e. the response 
of a single 2L atom, with respect to the Dopplcr shift kv 
is presented in Figs. [3] (b-d), red solid line, for different 
values of the area of D and EM pulses at a constant tem- 
poral duration of the pulses. It is clearly shown in the 
figure that, in the velocity domain pointed out above, 
with kv ~ (Ao + A) , the 2L atoms arc strongly per- 
turbed by the pulses and they do not end at the initial 
state W = -1. 

Therefore, the macroscopic response of the inhomo- 
geneously broadened 2L medium considered is provided 
predominantly by atoms in the velocity domain of nona- 
diabatic response. As known the macroscopic absorp- 
tion of atoms within a particular velocity domain of the 
Dopplcr profile is defined by absorption of individual 
atoms [represented by Wfi n , see red solid lines in Figs. [3] 
(b-d) and Fig. [5] (e) ] weighted by the probability den- 
sity for these atoms to belong to this particular velocity 
domain within the Maxwell distribution (shown by the 
dashed lines in Fig. [3] (a) and Fig. [U (e) ). Therefore, 
the absorption of the 2L medium considered could be re- 
duced either by decreasing the absorption of individual 
2L atoms or by moving the velocity domain of nonadi- 
abatic response to the wings of the Maxwell distribu- 
tion where the probability density for atomic velocities is 
strongly reduced. We found that by increasing the am- 
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FIG. 5: (Color online) Time evolution of the macroscopic 
population difference (W) and the medium coherence (U, V)) 
(b) in the inhomogeneously broadened two-level system un- 
der excitation of EM and D pulses shown in (a), (c) shows 
the evolution of the mixing angle, 9, for atoms with zero pro- 
jection of the velocity onto the light propagation direction. 
The pulse areas are An — Aa = 107r . Distribution of the 
final amplitude for the single atom variables U and W on the 
Doppler shift kv is shown in (d) and (e) , respectively. Dashed 
lines in (d) and (e) are plotted to visualize the Maxwell distri- 
bution of the atomic velocities. The maxima (upper line) and 
minima (lower line) for the variable 9 are shown in (f), where 
dots represent its zero final value. The other parameters are 
the same as in Fig. [3] 



plitudc A of the D pulse (increasing simultaneously the 
amplitude of the EM pulse to keep the same pulse areas) 
the domain of nonadiabatic response could be shifted to 
the wing of the Doppler profile, as it is clearly shown in 
Figs.[2Ib-d). In this way its negative role is reduced due 
to the much lower relative atomic population at the wing 
of the Maxwell distribution shown by the dashed line in 
Fig-E^a), almost all the atoms undergo an adiabatic pro- 
cess and, consequently, the medium becomes essentially 
transparent. 

In order to verify the previous statement, we present in 
Fig.[5]further results corresponding to the case considered 
in Fig.[3]^d), which describes an inhomogeneously broad- 
ened system of 2L atoms excited by D and EM pulses of 
the form given in with area of A = 107r. Fig. 02a) 
shows the sequence of the exciting D and EM pulses in 
such case. The time evolution of the macroscopic vari- 
ables (U, V, W) (obtained at each time step by the sum- 
mation of single atom response weighted in accordance 
with the Maxwell distribution of the atom velocity) and 
of the mixing angle 9 (the latter is presented only for 
atoms with kv — 0) for the ASIT process is plotted in 
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Fig. EJb) an d Fig. E{c), respectively. One can see that 
the ASIT phenomenon has predominantly a dispersive 
nature, since the quantity V takes values much smaller 
than those taken by U [Fig.[SJb)]. In Fig.EJthe amplitude 
of the D pulses has been chosen two times larger than the 
Dopplcr width of the inhomogcneous transition. In such 
a case, the spectral components of the EM pulse induce 
almost identical evolutions in 2L atoms from different ve- 
locity classes of the Maxwell distribution (except for the 
perturbed domain mentioned above at kv ~ (Ao + A) ). 
Indeed, we found that a given exciting pulse in the inho- 
mogeneously broadened 2L system will produce almost 
the same response that would be induced in a homoge- 
neously broadened 2L system. It should be mentioned 
that for small values of the amplitude of the D pulse, the 
perturbed domain is moved to the center of the Maxwell 
distribution, as shown in Fig. [3] This strongly modi- 
fies the response of the inhomogeneously broadened 2L 
atoms because large part of them are transferred to the 
excited state, making the ASIT process ineffective. 
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IV. ASIT VERSUS SIT 

In this section we compare the robustness of the ASIT 
technique with respect to the SIT one by studying the ef- 
fect of the variation of the pulse parameters. We consider 
the evolution of an inhomogeneously (Doppler) broad- 
ened 2L system. Fig. El shows the numerical results 
for the final values of the macroscopic population differ- 
ence Wfi n [Figs. El (a) and [6] (b)] and for the time aver- 
aged quantity g 2 ((QV) 2 ) / (f2 2 ) describing the absorp- 
tive character of the process [Fig. Etc) and Fig. Et[d)] , in 
the parameter plane defined by the area of the D and/or 
EM pulses and the pulse width. Fig. EI a ) and Fig. [6fc) 
correspond to a SIT process and Fig. [6](b) and Fig. Eld) 
correspond to an ASIT process. 

Since the ASIT process is based on a doublc-STIRAP 
technique, and since the STIRAP mechanism is more ef- 
fective when the exciting pulses have equal amplitudes, 
we have assumed that D and EM pulses have equal ar- 
eas, Aq = A&, and time durations in the simulations 
of Fig. El Moreover, it has been assumed that the delay 
time between consecutive D and EM pulses is At = 1.5r, 
for which it is known that STIRAP transfer by Gaussian 
pulses is most efficient . For the SIT process we have 
considered only an EM pulse of area Aq (Aa = 0). 

White regions in Fig. EIJa) and Fig. Eljb) correspond to 
the domains with Wfi„ — > — 1 , which means that the final 
state is the same as the initial one and thus there is trans- 
parency. Fig. EI a ) shows that to achieve transparency in 
the SIT case the pulse area needs to be precisely con- 
troled (multiples of 2tt), especially for short pulses. But 
even for these particular values of the pulse area, Doppler 
broadening and relaxation processes prevent from reach- 
ing complete transparency (the best values in Fig. Eta) 
are obtained for pulse area equal to 2n and short pulse 
duration, where Wf in ~ -0.973 at r n = 2 x I0 _3 7 _1 ). 



FIG. 6: Contour plot of the final values of the macroscopic 
population inversion, Wfin, after SIT (a) and ASIT (b) pro- 
cesses in a Doppler broadened 2L system as a function of the 
pulse area and duration. A a = An in (b) and (d), and A a ~ 
in (a) and (c). Other parameters of the 2L system and of the 
pulses are the same as in Fig. [3] White regions in (a) and (b) 
correspond to the domains where Wfi n — > — 1. Plots for time 
averaged values of g 2 ({OV) 2 ^ / (^ 2 ) in g 2 units indicate the 
pure absorptive nature of the SIT phenomenon (c) and the 
more dominantly dispersive nature of the ASIT process in the 
domain close to the left vertical axis (d). 



For longer pulses (and thus with smaller peak amplitude) 
the broadening due to Doppler effect (and also due to 
spontaneous emission) is more present because its rela- 
tive influence on the value of the generalized Rabi fre- 
quency increases. In contrast, Fig. E(b) shows that in 
the ASIT case transparency is feasible with pulses of ar- 
bitrary area, provided this area is larger than a value 
(approximately 47r) which is almost equal to that re- 
quired by the adiabaticity condition (J3]) in the doublc- 
STIRAP process. Note also in Fig. El^b) that the degree 
of transparency is larger for smaller pulse durations (i.e., 
near the vertical axis), where values very close to — 1 
{Wfi n ~ —0.995 at tq = 2 x 10 _3 7 _1 for pulses of area 
An = A/\ = 10-7T ) are achieved. For larger pulse dura- 
tions (and thus lower pulse amplitudes) it is progressively 
more difficult to push the velocity classes not fulfilling 
the adiabaticity condition (i.e., classes kv ~ (Ao + A)) 
to the wings of the Maxwell distribution and larger pulse 
areas are required for that. Also the influence of spon- 
taneous emission is larger for longer pulses. Note also 
in Fig. El^b) that near the horizontal axis there is a sec- 
ond region with values of Wfi„ close to -1. Such region 
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FIG. 7: Contour plot of macroscopic variable Wfin un- 
der conditions of ASIT in the plane (An = 4a, Awd) for 
a Doppler broadened 2L system of 87 Rb atoms excited by 
EM and D pulses of equal areas for two values of the 
time duration rn = ta = 0.0l7 _1 (a) and rn = ta = 
0.057 -1 (b). The considered transition of the 87 Rb atoms 
is (\5 2 S 1/2 ,F = 2,m F = 2) <-> |5 2 P 3/2 ,f = 3,m F = 3) with 
A = 780 nm, 7 = 5.83 x 2tt x 10 6 s' 1 and T = 27. 

does not correspond to pure ASIT, since the adiabaticity 
condition is not fulfilled; rather, it can be seen as a remi- 
niscence of the SIT effect due to the EM pulse (compare 
with Fig. Efa)), combined with the relaxation processes. 
Fig. [UJc) and Fig. [f^d), confirm the previous discussion. 
The absorptive character is smaller in the ASIT case than 
in the SIT case. 

It is known that the Doppler width in optical transi- 
tions can be modified by cooling or heating the atoms. 
In such cases it is possible to change the Doppler width 
in a large range of values from below to few orders above 
the homogeneous width. In order to show the feasi- 
bility of ASIT for real hot atoms we have investigated 
the effects of the Doppler width in specific atomic tran- 
sitions. The results are summarized in Fig. where 
we plot the dependence of Wfi„ in a hot gas of 87 Rb 
atoms as a function of the amplitudes of the EM and 
D pulses and of the Doppler width of the inhomoge- 
neously broadened transition (\5 2 Si/ 2 , F = 2,mp = 2) 
<-> \5 2 P 3 / 2 ,F = 3, nip = 3), for two values of the pulse 
width r = T A = O.OI7- 1 (Fig. Ha)) and = 0.057" 1 
(Fig. EJb)). It is clear from Fig. [7] that nonadiabatic 
effects arc negligible for pulses with large area and short 
time duration. In contrast they are pronounced for long 
pulses and also for short pulses with small area, when the 
amplitude A of the D pulse becomes comparable to the 



Doppler width of the transition. By taking into account 
that at room temperature (T = 300°K), the Doppler 
width is equal to Au>d — 5307 we can conclude, [see Fig. 
[7{a)] , that the reported phenomenon of ASIT can be fea- 
sible experimentally in an atomic vapour cell of rubidium 
at room temperature and under excitation of nanosecond 
EM and D pulses with pulse area of few it. 

In conclusion, we have demonstrated the feasibility 
of the adiabatic self induced transparency (ASIT) tech- 
nique, which yields Doppler-free transparency for electro- 
magnetic pulses in a 2L system. Our technique consists of 
driving the 2L system with a laser pulse whose detuning 
with respect to the atomic transition and whose EM field 
amplitude are temporally varied in a proper shape. The 
detuning modulation can be introduced cither by using 
the Stark or the Zceman effect to obtain timc-dcpcndcnt 
energy levels in the atom/molecule, which gives rise to 
time-dependent detunings, or alternatively by introduc- 
ing a suitable chirp on the EM pulse. We have studied the 
robustness of the reported phenomena against modifica- 
tions of the pulse parameters and we have shown that the 
ASIT process is much more robust than the conventional 
SIT, in the sense that the area of the pulses needs not to 
be adjusted with precision and it can be much less sensi- 
tive to Doppler broadening. Because of the adiabaticity 
condition, the pulse areas for ASIT need in general to 
be somewhat larger than for SIT, for which a 2n area 
is sufficient, but it could be implemented with no essen- 
tial difficulty in real Doppler-broadened atomic systems 
such as for instance rubidium hot vapors. Although the 
results presented in this article have been obtained for 
Gaussian pulses we have also checked the validity of the 
ASIT technique for pulses with other forms of temporal 
profiles. 



Acknowledgments 

We acknowledge support from the Spanish Ministry 
of Education and Science (MEC) under the Programa 
Ramon y Cajal and under the contracts FIS2005-07931- 
C03-03, FIS2005-01497 and FIS2008-02425, from the 
project "Quantum Optical Information Technologies" 
within the Consolidcr-Ingcnio 2010 program under the 
contract CSD2006-00019. Support from the Catalan 
Government under Ci^-ED-program and under the con- 
tracts SGR2005-00358 and SGR2005-00457 is also ac- 
knowledged. 



[1] A.Kasapi et ai, Phys. Rev. Lett. 74, 2447 (1995); L.V. 

Hau et al, Nature (London) 397, 594 (1999). 
[2] M. Xiao et ai, Phys. Rev. Lett. 74, 666 (1995); 

O.Schmidt et ai, Phys. Rev. A 53, R27 (1996); 

M.M.Kash et ai, Phys. Rev. Lett. 82, 5229 (1999). 



[3] Fleischhauer and M.D. Lukin, Phys. Rev. Lett. 84, 5094 
(2000); M.D. Lukin, S.F. Yelin, and M. Fleischhauer, 
Phys. Rev. Lett. 84, 4232 (2000); M. Fleischhauer and 
M.D. Lukin, Phys. Rev. A 65, 022314 (2002). 

[4] J.I.Cirac, et ai, Phys. Rev. Lett. 78, 3221 (1997); C. Liu 



8 



et al, Nature (London) 409, 490 (2001); D.F. Phillips et 
al, Phys. Rev. Lett. 86, 783 (2001); M. D. Lukin, Rev. 
Mod. Phys. 75, 457 (2003); A.Kuzmich, et al, Nature 
(London) 423, 731 (2003); C.H.van der Wal, et al, Sci- 
ence 301, 196 (2003); M.D. Eisaman et al, Nature 438, 
837-841 (2005). 

[5] A.V. Turukhin et al, Phys. Rev. Lett. 88, 023602 (2001); 
J. J. Longdell et al, Phys. Rev. Lett. 95, 063601 (2005). 

[6] J. Q. Liang et al, Phys. Rev. A 65, 031801 (2002); M. S. 
Bigelow et al, Phys. Rev. Lett. 90, 113903 (2003); M. S. 
Bigelow et al, Science 301, 200 (2003); P.-Ch. Ku et al, 
Opt. Lett. 29, 2291 (2004); E. Baldit et al, Phys. Rev. 
Lett. 95, 143601 (2005). 

[7] M. Scalora et al, Phys. Rev. E 54, R1078 (1996); 
M.Notomi et al, Phys. Rev. Lett. 87, 253902 (2001); 
H. Gersen et al, Phys. Rev. Lett. 94, 073903 (2005); 
Yu.A.Vlasov, Nature (London) 438, 65 (2005). 

[8] D.J.Joannopoulos, R.D.Meade, J.N.Winn, Photonic 
Crystals (Princeton University Press, 1995). B.E.A.Saleh 
and M.C. Teich, Fundamentals of Photonics (Wiley & 
Sons, Ney York Chichested, 1991). R.E.Slusher, and 
B.J.Eggleton (eds.) Springer series in Photonics: Non- 
linear Photonic Crystals (Springer, Ney York, 2003). 

[9] K.Y. Song et al, Opt. Lett. 30, 1782 (2005); Y.Okawachi 



et al, Phys. Rev. Lett. 94, 153902 (2005); J.E. Sharping 

et al, Optics Exppres 13, 6092 (2005). 
[10] F. Silva, et al, Europhys. Lett. 51 (3), 286 (2000); F. 

Silva, et al, Phys. Rev. A 64, 033802 (2001). 
[11] E. Kuznetsova, O. Kocharovskaya, P. Hemmer, and M. 

O. Scully, Phys. Rev. A 66 063802 (2002). A. L. Alexan- 
der et a/.. larXiv:0806.4324l 
[12] S. L. McCall, E. L. Hahn, Phys. Rev. 183, 457 (1969); 

G. L. Lamb, JR., Rev. Mod. Phys. 43, 99 (1971). 
[13] E.Arimondo, in Progress in Optics XXXV (ed. EWolf), 

259 (Elsevier Science, Amsterdam, 1996). 
[14] S.E. Harris, Phys. Today 50, 36 (1997); M.O. Scully and 

M.S. Zubairy, Quantum Optics (Cambridge University 

Press, Cambridge, U.K., 1997); J. P. Marangos, J. Mod. 

Opt. 45, 471 (1998); M. Fleischhauer, A. Imamoglu, and 

J. P. Marangos, Rev. Mod. Phys. 77, 633 (2005). 
[15] K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. 

Phys. 70, 1003 (1998). 
[16] N. V. Vitanov, et al, Annu. Rev. Phys. Chem. 52, 763 

(2001). 

[17] N. V. Vitanov, and B. W. Shore, Phys. Rev. A 73, 053402 
(2006). 



